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ABSTRACT 

A number of recent surveys for gravitational lenses have found examples of double Einstein 
rings. Here, we investigate analytically the occurence of multiple Einstein rings. We prove, 
under very general assumptions, that at most one Einstein ring can arise from a mass distri- 
bution in a single plane lensing a single background source. Two or more Einstein rings can 
therefore only occur in multi-plane lensing. Surprisingly, we show that it is possible for a 
single source to produce more than one Einstein ring. If two point masses, or two isothermal 
spheres, in different planes are aligned with observer and source on the optical axis, we show 
that there are up to three Einstein rings. We also discuss the image morphologies for these 
two models if axisymmetry is broken, and give the first instances of magnification invariants 
in the case of two lens planes. 

Key words: Gravitational lensing 
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1 INTRODUCTION 

In his seminal article on the gravitational lensing effect, Einstein 
(1936) discussed the circular image of a point-like source, noting 
that "there is no hope of observing this phenomenon directly." But, 
thanks to advances of instrumentation since then, arcs of partial 
and even complete Einstein rings are now being routinely found. 
Recently, surveys have found the first instances of multiple Ein- 
stein rings. This includes the partial double Einstein rings of SDSS 
10924 +0219 discovered from HST images by the COSMOGRAIL 
team l lEigenbrod et al. 20061), and of SDSS J094 6+1006 found by 
the Sloan Lens ACS Survey dOavazzi et al. 20081) . The Cambridge 
Sloan Survey of Wide Arcs in the Sky (CASSOWARY, Belokurov 
et al. 2008) has also uncovered a number of examples of multiple 
ring systems, such as CASSOWARY 2. Such systems may offer 
valuable insights into the mass distribution of the lensing galaxies. 

From a theoretical point of view, the possibility of forming 
multiple Einstein rings has long been known in the case of the 
strong deflection limit near the photon sphere of a black hole (for a 
review, see e.g. Nemiroff (1993) and references therein). However, 
a systematic investigation for the weak deflection limit appears to 
be lacking so far, and we present results to this end for the case of 
one and two lens planes in this paper. After an outline of the gen- 
eral lensing setup, the condition for Einstein rings is derived in §2. 
In the single lens plane case considered in §3, we prove that, under 
rather general assumptions, multiple Einstein rings cannot arise. 
We therefore proceed to two lens planes in §4 and consider two 
singular isothermal spheres and two point lenses as simple models 
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where multiple Einstein rings of a single source do, in fact, occur. 
In the first example, there are up to two Einstein rings due to the 
singular isothermal spheres in different planes, and another ring if 
the second lens is also luminous. Similarly, two point lenses can 
give rise to up to three Einstein rings overall. Therefore, the usual 
supposition that arcs of multiple Einstein rings indicate the pres- 
ence of as many sources at different distances is not necessarily 
correct. We briefly discuss the image configuration in these cases if 
axisymmetry is broken. For models with two lens planes, we also 
find that analogues of the invariants of the signed magnification 
sum (see e.g., Witt & Mao 1995; Hunter & Evans 2001) hold in the 
domains of maximal image multiplicity. 

Regarding notation, we write and A^, for the gradient 
operator and Laplacian, respectively, expressed in the same co- 
ordinate system as the vector v. Furthermore, is the vec- 
tor norm with respect to the Euclidean metric, and square brack- 
ets [u] denote functional dependence on the variable u. The set 
M.Q = {a; G R : 2: ^ 0} denotes the set of all non negative re- 
als. Subscripts can label vector components or images according to 
context. The universal gravitation constant and the speed of light in 
vacuum are represented by G and c, respectively, as usual. 



2 LENSING FRAMEWORK 
2.1 General setup 

Gravitational lensing in the weak deflection limit is conveniently 
described in terms of the impulse approximation, with piecewise 
straight light rays in flat space between the source, lens planes 
and observer (for a comprehensive introduction see e.g. Schnei- 
der, Ehlers & Falco 1999). We consider a point source in the source 
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plane S and two lens planes -L'^\ -L'^' with angular diameter dis- 
tances D\ ^ D2 ^ Ds between the observer and L'^^\ L'^', S, 
also Di2,Dis between L'^^ and L'^' , 5", and D2s between L*^^ , S, 
respectively. This setup is illustrated in figure [T] Using Cartesian 
coordinates 77 in 5 and $y\ $}'^'^ in 1/'^', L'^^ measured from 
some optical axis, and for the deflection angles of the 

light rays crossing the respective lens planes, the lens equations be- 
come. 



(2)rd:(2)l 



(1) 

(2) 



Assuming that the optical axis passes through the centre of the 
projected mass distribution in L'^', the deflection potentials in 
jr^(i) ^ £^(2) g^jj written with an overall scaling factor proportional 
to the lens mass such that 

^(i)[^(i)] = m(^V*''K<''], (3) 

q>(2)[^(2)] ^ M(2)/(2)[^(2)_^(2)]^ (4) 

where the centre of the projected mass distribution in L'^^ may 
be taken as = (^^^' , 0) without loss of generality. Then the 
deflection angles are 



(0 rtf(«)i 



(5) 



In this notation, the deflection potentials are solutions of Pois 
son's equation, corresponding to Einstein's equation in this quasi 
Newtonian approximation, 

S-kG 

„2 



-e«kW], ze{i,2}, 



(6) 



where E'"' is the projected surface mass density in the ith lens 
plane. Given a model for the surface density E'"', the integration 
constants for the corresponding deflection potential 5'''' are chosen 
such that the norm of the deflection angle at infinity is as small as 
possible, that is, zero for realistic lens models. 

While ^'^^ / Di would correspond to the usual angular coordi- 



nate in I/' we shall find the space of parameters more conven ient 
if the follow 

i2) e>/D2 „ rj/D. 



.(1) 



(2) ^ e'^/D2 



(7) 



with parameters 

^(1) +^{2)^1/2 



(8) 



where ^i^^) = Af^^^ and ^t^^^ = M'~'^\D2s / D2) and 

D2sD^ 



D^sD2 



DisD2 



(9) 



Here, < /? < 1. In particular, P — if and only if D12 = and 
/3 = 1 if and only if D2S = 0, provided that Di / 0, D2 / and 
Ds 7^ 0. Then, the lens equations ([TJ, ([2} can be rewritten thus 



.(2) 



,(1) 
.(1) 



(1)/ 
(1)/ 



introducing the normalized mass parameters 



,(1) 



,(1) 
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Figure 1. Multiple Einstein rings. The source plane S and the two lens 
planes L^^\ L'^' are at angulai' diameter distances Ds, Di, D2 from the 
observer O, respectively. A point source in S and a luminous lens in L^^^ 
on the optical axis (solid line) a can give rise to three Einstein rings with 
angular radii 8^; 1,9^ 2^^E 3 as discussed in sections 4.1.2 and 4.2.2. 
Schematic light rays ai'e shown as dashed hnes. 



.(2) 



,(2) 



(13) 



such that m<^^ + m^^^ 



2.2 Einstein rings 



Einstein rings, within the geometrical optics approximation of the 
standard lensing framework, are infinitely magnified, circularly 

symmetric images. We therefore stipulate circular symmetry about 

f 2) 

the optical axis such that x)^{ =Q and 



= n 



G {1,2}, and \\x 



is used for notational simplicity. But then the lens equation jlOt 
implies that a;'^-* = x*-^^ [^^'^'l and i ll It becomes 



y = 



,(1). 



(1) 



m 



xW dxW 



1 - I3m^ 



df 



(2) 



.(1)1 



Introducing plane polar coordinates {x^^\( 
minant of the lensing map ( I14t is therefore 



(14) 

, the Jacobian deter- 



det J 



1 

7w 



Syi 

axW 

9y2 



dy2 



F { F + x 



(1) 



dF 



(15) 



^ ^ ^ dxW 

The condition for critical curves det J = gives rise to two classes 
of critical curves. The first, given by f [a;''^-'] = 0, defines tangen- 
tial critical circles in L'^^ which are solutions of the lens equation 
( 114b , and hence infinitely magnified images, and map to the caustic 
point y = 0. The other solution F+x^-^^dF/dx^^^ — gives radial 
critical circles in L^^^ which map to caustic c ircles and hence define 
domains of constant image multiplicity in S jSchneider et al. 19991 
p. 233). 

Definition 2.1. An Einstein ring is a circular, critical image of a 
point source at y = whose radius is a solution of f [a;''^-'] = 
0, a;(i) > 0. 

Finally, we note that if y 7^ so that the axisymmetry is broken, 
discrete images are forme d at some a;*^^^ which have finite signed 
magnification [Schneider et al. I999I p. 162) 
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detJ[x(i)]' 



(16) 



3 ONE LENS PLANE 

We continue by specializing the previous discussion to the simpler 
case of a single lens plane L*-^^ — L^'^^ such that m^^-* = 0, /3 = 0. 
Notice, then, that l[5), l|7J and l[8j yield 



DiD 



(1) ^ DiDi.M^^^ 



1« A (1) 



which is the normali zed deflection angle in the standard form 
jSchneider et al. 19991 p. 158). Similarly, the Poisson equation ^ 
becomes 



(1) 



BttG, 



(1) 



2k 



with the usual definition of the critical surface density Ecrit 



c Di/dnGDi Di 9). Using definition l2.1l and writing ||q;|| = a = 
df^^^/dx^^\ the problem of finding Einstein rings therefore re- 
duces to a fixed point equation 



(1)1 



.(1) ^(1) 



>0, 



(17) 



subject to the Poisson equation for the given mass distribution of 
the lens, 



(18) 



Hence, we need to define general, yet astrophysically sensible lens 
models k[x-'^' ] that allow for Einstein rings. Apart from the circular 
symmetry inherent in the problem, it would be plausible to stipulate 
that K decreases monotonically with x'^^K However, we shall use a 
condition even weaker than monotonicity, namely that, at every ra- 
dius, the surface density be smaller than the average density of the 
mass enclosed. This is natural for self-gravitating and hence cen- 
trally condensed systems, and a more specific model for the mass 
distribution need not be assumed. 

Definition 3.1. The gravitational lens is defined by a normalized 
surface density k as mass model such that the following conditions 
are fulfilled. 

(i) Continuity: k : L^^^ be a continuous function except 
at e for singular lenses. 

(ii) Circular symmetry: ft = only. 

(iii) Finiteness: k < 00 for x'^^ > 0, fi:[0] = 1/Ci with constant 
Ci 0, and liiiij;(i)_,oo ^2;'^^ = C2 with constant ^ C2 < 00. 

(iv) Self-gravitation: /tfa;^'"-'] < /t[a::'"'-'] where ^[2:'"'"^] = 

(^(i))2 /o K\x\xdx. 

For example, the point lens has Ci =0, C2 = 0, and the isother- 
mal sphere Ci = 0, C2 > in our notation. Both are usually 
regarded as singular, since their projected surface densities diverge 
at the centre. Moreover, the total mass is infinite in the isothermal 
case, and we should expect C2 = for more realistic lenses. How- 
ever, it turns out that we do not need to require this for our purposes 
and hence define singular lenses simply as follows. 

Definition 3.2. A gravitational lens is called singular if, and only 
if, Ci = 0, and is called non-singular otherwise. 

With these definitions, one can now prove that there is at most one 
Einstein ring. This result follows from existence, whose necessary 



and sufficient condition is established in theorem [3.1l and unique- 
ness, shown in theorem [3 . 2 1 below. First of all, however, we shall 
need a lemma regarding the surface density, and a lemma concern- 
ing the deflection angle. 

Lemma 3.1. Given a gravitational lens in the sense of definitions 

(1)1 ^ 



[33] and [221 then the normalized surface density fulfils k[x 

«:[0] Va;(i> > 0. 

Proof. If the lens is singular, then the lemma follows immediately 
from definition l3.ir iii). Otherwise, by definition l3.ir iv), we have 

«:[a;('>] < k[x^^^] < max ac[x] Vx*''. 

Suppose that this maximum is attained at some x' < x^^\ Then 
the previous inequalities hold for k[x'] as well, with maximum at 
x" < x', say. Repeating this argument sufficiently often shows that 
k[x'^-'] < k[0], as required. □ 

Lemma 3.2. Given a gravitational lens in the sense of definitions 
l3.1l and lT2] then the normalized deflection angle a has the follow- 
ing properties. 

(i) Smoothness: a : i'^' M.q is a smooth function except at 
e L*^' for singular lenses. 

(ii) Circular symmetry: a — afx'^-*] only. 

(iii) Properties at the centre: For non-singular lenses, a[0] = 
and the derivative da/da; ' ^ ■* [0] = k[0], for singular lenses a[0\ > 
0. 

(iv) Finiteness: lim (i)^ a — A with constant A < 00. 



Proof. Definition [2211) implies property (i), and property (ii) fol- 
lows immediately from definition (STTJii) and the uniqueness of so- 
lutions of Poisson's equation. Circular symmetry and smoothness 
at the centre imply that a[0] — for non-singular lenses, and a[0] 
is some positive, possibly infinite, value for singular lenses. Now 
integrate Poisson's equation 1 II8I 1 to find 

l^(a:r(i)) ^ (19) 



2k 



.{1)1 



x(i) dxW 
2 ..(1) 



x]xdx. 



(20) 



using a[0] = 0. Then l |19l l implies the second property of non- 
singular lenses in (iii), 

r-(l)^2,\ 



lim 



2 ^(1)^0 \ dxW 
da 



da 



r(i) 



using again a[0] — 0, and property (iv) follows from equation jlOl 



and definition l3.U iii) since 
lim a — 2 lim kx'^^ = 2C2, 

^(1) — voo a;(^) — ^00 

SO y4 = 2C2 < 00, as required. This is also true for singular lenses 
where the integration of il9l has to start at some x > 0, so that 
the limit for a as x^^^ 00 is the same as before plus some finite 
integration constant. □ 

Theorem 3.1. (Existence) Given a gravitational lens in the sense 
of definitions 13. 1 1 and 13.21 then at least one Einstein ring exists if, 
and only if, the condition k[0] > 1 holds. 

Proof. First, we prove the sufficiency of the condition: k[0] > 1 
3 Einstein ring. Let us consider the graph of the deflection angle 
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graph a = {x^'-> , a(x^'->)) e R,j x KJ. 

Then, according to equation l ll7t . the existence of an Einstein ring 
is equivalent to the existence of a fixed point and hence an inter- 
section of graph a with the diagonal (x'^-*, a;'^-*) G 
some s'^^ > 0. Let Ai = {{x,y) £ 



at 



y > x'^^} and 

A2 = {{x,y) e X : y < a;'^'} be two domains to the 
left and right of the diagonal, respectively. For non-singular lenses, 
lemma [X2l iii) implies that a;^^-* = is a fixed point but no Einstein 
ring, and k[0] > 1 da / dx'-^^ [0] > 1 graph a £ Ai for 
s'^^ — > 0. This is also true for singular lenses because here a[0] > 
0. On the other hand, lemma [X2l iv) implies that graph a £ A2 for 
x-'^' — > 00. This ensures the existence of at least one fixed point 
with x'^' > and hence Einstein ring. 

Now we show necessity: 3 Einstein ring k[0] > 1. This 
statement is equivalent to its contraposition, k[0] ^ 1 ^ $ Ein- 
stein ring, which we prove by contradiction. So suppose k[0] ^5 1 



and 3 Einstein ring, then there is some x^' > such that a[a;'^'] — 
a;^' by equation l |17| l. By recasting Poisson's equation l[T8) and us- 



.(1)1 



ing lemma 1371 



.(1)1 



2k[x 



(1)1 

E 1 



1 < 2k[( 



< 



da 
da 

since by assumption k[0] ^ 1. But on the other hand, for a fixed 
point a;^' > to exist, we must require dQ/da::'^-'[a;^^] ^ 1. 
To see this, notice first of all that the assumption «:[0] ^ 1 also 
means, by definitions IS.U iii) and 13.21 that we only need to con- 
sider non-singular lenses here. Now, by lemma [X2l iii). this implies 
that da/dx'^) [0] s: 1. By differentiating ^9^ and using THopital's 
rule. 



d^a 



[0] 



4 dK 



< by lemma [37T 



dx(i)2^ 3dx(i) ' 

so even if k[0] = 1, graph a £ A2 for small x^^\ and the result 
follows. This completes the proof by contradiction and hence the 
proof of the theorem. □ 

Theorem 3.2. (Uniqueness) Given a gravitational lens in the 
sense of definitions |3.1| and |3.2| then if an Einstein ring exists, there 
is exactly one. 

Proof. The existence of an Einstein ring E implies that k[0] > 1 
by theorem [3T| and that there is some x^' > such that a[x^''] — 
by equation illl . Using the integral equation l l20b and defini- 
tionlTTliv), we obtain 



a[xW] 



.(1) 



(xW)2 Jo 



K[x]xdx = fi:[x^''] 



Hence by the Poisson equation il9\ and definition l3.ir iv). 

da 
d^ 
da 



= 2k[x'-^'>] - 1 < 2R.[x^^'>] 



dx(i) J 



< 



there is another Einstein ring E' at x^V > a;^'. But by the 
same token used in the sufficiency proof of the previous theorem, 
graph a £ Ai for small x'^-* and graph a £ A2 for x*-^-* — > 00, 
so we need da/dx'^''[x^,^] ^ 1 for the fixed point of E' to exist. 
This cannot be according to l l21b . and the result follows. □ 



Since multiple Einstein rings cannot occur in a one lens plane set- 
ting by theorems l3.1l and [T2l we now discuss to two simple models 
in the two lens plane case. 



4 TWO LENS PLANES 

4.1 Singular isothermal spheres 

4.1.1 Lens equation 

The first system of two lenses discussed here consists of a singu- 
lar isothermal sphere both in L^^' and L'^^\ This is hence a spe- 
cial case of the cored, spherically symmetric lenses in two lens 
planes considered by Kochanek & Apostolakis (1988) in the con- 
text of a numerical study of lensing cross-sections. However, one 
can easily extend this model, which has two identical lenses, to 
include the effect of different surface densities. Because of the cir- 
cular symmetry required for Einstein rings, we shall assume the 
two singular isothermal spheres to be centred on the optical axis 
so that = aga in, and consider projected surface densities 



jSchneider et al. 19991 p. 243) 



2GC(') ' 



G{1,2}, 



where = H^'''!! and (j^''-* is line of sight velocity dispersion 
of the isothermal sphere in the ith lens plane. Now according to 
Poisson's equation the deflection potentials l|3j, & become 



^M^^'i^'' where Af 
and hence the lens equation l ll4t . 



47r(j 



2(i) 



ie{l,2}, 



(1) /, rn± 
V xl 



J) 



(21) 



in which m'^' , mP'' from equations ( I12t and l ll3t have been com- 
bined to define a new mass parameter in this case. 



m± 



(1) I D2 (2) 

1^ ' ± -:—m} ' 



(22) 



where the positive sign is valid for x'^' P^om^^^ and the nega- 
tive sign for x'"'' < f3^o'm''-^K 



4.1.2 Einstein rings 

Now if y = 0, F[x^-^''] = 1 — m±/x^^\ no radial critical circles 
exist. Definition 12.11 and | |2U imply that the radii of Einstein rings 
are given by 



(1) 

X — m±, 



(23) 



subject to the two domains of the mass parameter l l22t . It turns out, 
then, that the point source in S lensed by both isothermal spheres 
in L'^' and L'^' always produces one Einstein ring of radius 



.(1) 



1 M^^>£)i, +M^^>_D2 



(24) 



and, in addition, provided M^^^/Af'^' < D2/D1, another one of 
radius 



^(1) _ 

■''E,2 — 



1 M'-^'> Dis - M'-'^^ D2 



(25) 
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Furthermore, there is a third Einstein ring of the isothermal sphere 
in L^^-* lensed by the one in L^^\ assuming, of course, that at least 
the former is luminous. Because the lens in L'^^ is extended, this 
second Einstein ring is not a circle but has some radial width. We 
therefore define the radius of this third Einstein ring x^^\ in 
to be that of the circular image produced by the centre of the lens 
(the most luminous part) in _L'^' . Its radius can be read off immedi- 
ately from the previous results by letting Ad^^^ = and replacing 
Ds, Dis by D2, D12, respectively, so that there are two different 
definitions of one for the two rings due to the point source, and 
one for the ring due to the isothermal sphere in Z/'^\ Hence, using 
llSJ and l|22j. 



D2 



1/2 



In view of observational applications, it is more convenient to ex- 
press these Einstein ring radii in terms of angular coordinates in 
L^^' . Given the small angles approximation, one can take the an- 
gular radius to be ^ = ^'^^Di = x'^'^o/^^i using Q. Hence for 
the three Einstein rings. 



M^^^ Dis - M^^^ D2s £2 



-M 



(1) 



(26) 
(27) 
(28) 



4.1.3 Image configuration 

According to the lens equation l l21b . this system of two singular 
isothermal spheres at different distances is a modification of the 
well-known single lens plane case in the sense that the two in- 
stances m± have to be distinguished here. It turns out, then, that 
up to four images can be obtained in the present case. To be more 
precise, consider the image configurations if axisymmetry is bro- 
ken such that y 7^ 0. Without loss of generality, one may take 
2/1 = y > 0, 2/2 = so that y = {y, 0). Images are hence collinear 
with the centre of i'^-* and the source such that x'^^ = (x-,0), 
say. Therefore \m and i22\ imply the following image positions 
Xi, i G {1, . . . , 4}, for given domains of y £ S, 



Xi 



y + 771+ > for 



y ^ -m+ + /3^om'^' = 



(29) 



Co D2s M'^'Di +M(')L>2 



Di D2 IJ.'-^^ + 
X2 = y + m->0 for (30) 

Co D2s M'^^'iDi- M'^^'>D2 



y < -m^ + l3^om^' = - 



Di D2 + 
y - m+ < for (31) 

(1) _ Co D2, M(i)Di +M(2)D2 



Di D2 + m''^ 

Xi — y — m_ < for (32) 

(1) Co D2s W-^^Di - Af(2>D2 



Dx D2 



The condition on y for M9\ to hold is clearly always fulfilled, and 
the condition for ( |31t defines a cut circle in S within which another 
image occurs. Hence the images given by x-i and x-i correspond to 
the single lens plane case if /5 = (see e.g. Schneider et al. 1999, 
p. 244). Similarly, the condition on y in l |30t defines a second cut 
circle. But because we need to keep y positive by setup, this cut 



circle only exists if M'^^ < M^^^D2/D 1, the same condition as 
for the second Einstein ring i27\ discussed in the previous section. 
Furthermore, if the source is outside of the second cut circle, then 
one image Xi now occurs on the opposite (i.e., negative) side of 
Z/'^' . If the source is inside the second cut circle, another image X2 
appears on the same side of L*-^' as the source. Finally, if j/ ^ 0, 
then all image positions approximate to the Einstein rings given 
by l l23t . as expected. However, it is interesting that, unlike in the 
single lens plane case, the cut circles and the Einstein rings do not 
coincide here. 

Now assume that y is in a maximal domain in S such that 
all four possible images Xi | |29I32| | are present. Then, given that 
detJ = l-m±/x<^' by (Bll, we can directly evaluate the sum of 
the signed magnifications l ll6t of the images to find 



(33) 



This is therefore a magnification invariant in the two lens plane 
case. 



4.2 Point lenses 

4.2.1 Lens equation 

We now turn to two point lenses. A thorough study of the caus- 
tic structure for two point lenses in general positions within three- 
dimensional space was done by Erdl & Schneider (1993). However, 
the aspect of multiple Einstein rings was not studied there explic- 
itly, so we present results to this end below. Using a setup analo- 
gous to the one in the previous section, consider two point masses 
Ml, M2 on the optical axis in L^^' , L^^ \ respectively. In this c ase, 
then, the surface densities are given by dSchneider et al. 199^ p. 
239) 

^ Ai^5'"-'^[^(')]^ i e {L2}, 

where i'^'*' denotes the two-dimensional delta function of the ith 
lens plane. Hence, the corresponding deflection potentials 



(i) 1 



(On 



Co 



where A/''' — 



4GM, 



follow from l|6j, and the lens equation l ll4t becomes 



y = x 



(1) 



1 - 



6 {1,2}, 



(34) 



(a;(i))2 -/3m(i), 

using ([Hi, lO and ^ jErdl & Schneider 1993h . Recall also that 
the range of the parameters used in ( I34t is ^ /3, m'^', m*-^' ^ 1 
by definition. 



4.2.2 Einstein rings 

Again, according to definition |2.1| Einstein rings are given by pos- 
itive solutions of F{x''^^) — for y — 0, which is a quartic in 



,{i) 



by lens equation i34l . Hence, letting (: 



.{l)^2 



r, we seek 



solutions 



F{r) 



,(1) 



,(2) 



1 



0, r > 0, 



r r — /3m(-'-) 
and obtain the squares of Einstein ring radii 



te = 



i ^1 + /3?n'') ± y^(l+/?m(i))2-4/3(m(i))2^ 



(35) 
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using l fT2t and l fT?t . Notice that 

(l + /3m('>)' > 4/3(m('>)' 
unless P — 1 A m'^^ = 1, and 



(36) 



unless /9 = V m 



(1) 



0. 



Therefore, in the general case < /3 < 1, < m'^^ < 1, there are 
exactly two Einstein rings with radii x^'j, a;^ j from l l35b due to a 
point source in S and the point lenses in L'^' and L^^-* . Otherwise, 
one ring becomes a critical point on the optical axis, and we recover 
the case of a single Einstein ring, as expected. 

If, in addition, the point lens in L'^^ is also luminous, it pro- 
duces a third Einstein ring at radius a-^ ^ due to the point lens 
in L^^K As before, this radius can be obtained simply by setting 
m''^^ = and identifying Ds , -Dis with D2 , -D12, respectively. The 
corresponding angular radii are again given 
x^^^^o/Di, and we find 



Oe,: 



Dis 
DiDs 



D12 



M(i), with mass parameter 



(37) 

(38) 
(39) 



using {Hi, lUl, (Hi and i fTSl l. 



4.2.3 Image configuration 



Now consider the case when y 7^ 0. As before, we can choose 
without loss of generality y = (j/,0), y > 0. The images are 
coUinear with the source so that x^^^ = (a-, 0). Hence solutions of 
the lens equation J34t are defined by a quintic in x. 



: 



yx 



(1 + I3m'-^'>)x^ + yl3m'-^'>x'^ + l3{m'-^'>fx. 



We immediately have a solution on the optical axis, X5 — 0, but 
this image is infinitely demagnified, and therefore invisible, since 
its signed magnification is given by 



MF5 



: because lim _F[(a::,0)] — —00 



det J[x5 

using l ll6t . The remaining roots Xi, i £ {1, . . . , 4}, are solutions 
of the quartic equation 

= a;"* - 1/2:^ - (1 + /3m*^')a;^ + j//3m'^'a: + /3(m'^')^ (40) 
We note that this immediately implies that the roots satisfy 



E 



Xi = y. 



Now Descartes' rule of signs shows that i40l has either two positive 
and two negative real roots, two positive real roots and a complex 
conjugate pair, two negative real roots and a complex conjugate 
pair, or two complex conjugate pairs of roots. So in the first case, 
there are four images, in the second and third case, two images 
each, and in the last case, none. As usual, these cases can be distin- 
guished by means of the cubic resolvent of i40l . With the standard 



substitution z = a; — ?//4 we can eliminate the cubic term in ( I40t 
to obtain a reduced quartic in z, 

= 2'* + Pz^ + Qz + R, where 

P - -^-(l + /3m(i)), 

Q = -^-f(l-/3m«), 

Hence the resolvent of l|40t is given by the cubic equation 
teronshtein & Semendvavev 1985L p. 121) 

= + Az'^ + Bz + C (41) 

where 

A^2P, B ^ -4R, C = -Q^. 
Notice, then, that A < and 

B = ^ + + (1 + /3m<^))' - 4/3{m^'^f > 
16 

using l l36t . Applying Descartes' rule of signs to | |4U with z re- 
placed by —2; shows that the cubic resolvent has no negative real 
roots, and so the quartic l |40| l does not have two complex conjugate 
pairs. Thus, there are always im ages present. Furthermore, one can 
consider the discriminant Dies dSronshtein & Semendvavev 1981 
p. 120) of the cubic resolvent to find 



- A 



2\ 3 



9 



/ 2A^ - 9AB + 27C 
V 54 



s; 0. 



Therefore the cubic resolvent has three positive real roots, and 
hence the quartic i40l has four real roots corresponding to four im- 
ages which are, in general, distinct. In the limiting case Dros = 0, 
however, double roots occur and we recover the single lens plane 
case with two images. This result is in agreement with Erdl & 
Schneider (1993) who have shown, using catastrophe theory, that 
two point lenses can produce either four or six images. This has 
also been proven by Fetters (1995) using Morse theory. It is clear, 
then, that the caustic domain giving rise to six images does not 
occur in our axisymmetrical case. The limiting behaviour of these 
four images can be read off directly from the quartic equation \40h 
For y ^ 0, there are two images with positive x and two images 
with negative x approaching the two Einstein circles. For y ^ 00 
the four images Xi satisfy 



a;i 



y, X2 



1,(1). 



xs 



,.(1). 



X4, 



Using l ll6t in these cases, we find that X2,X3,X4, become infinitely 
demagnified so that only xi with signed magnification jJ.[xi] — » 1 
remains as j/ — > oo, as expected. 

Furthermore, it turns out that the images obey an invariant of 
the signed magnification sum. To see this, we can re gard the presen t 
case as an example of a more general argument 1 Werner 2007h . 
First, notice that all five possible roots of the lens equation ( |34| > 
for y > are real by the previous discussion such that all of S\0 
is a maximal domain. Also, the deflection angle tends to zero as 
y ^ 00. Then the lens equation can be complexified such that im- 
ages correspond to fixed points of a complex rational lensing map. 
This in turn induces a map on complex projective space which 
is holomorphic almost everywhere and, in particular, at the fixed 
points. Then the holomorphic Lefschetz fixed point formula can be 
applied to yield 
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5 

Y,^l[x,] = l. (42) 

i = l 

In fact, this statement is also true for the four visible images be- 
cause xs is infinitely demagnified as noted above. 



5 CONCLUSIONS 

In this article, we have presented an analytical study of multiple 
Einstein rings in the weak deflection limit of strong lensing, and 
now summarize the three main results. First, it was proven gener- 
ally that at most one Einstein ring can occur for one lens plane. A 
natural and weak assumption was used here, namely that the nor- 
malized surface density always be smaller than the average surface 
density. 

Accordingly, we turned to models with lenses in two planes, 
and in which source, observer and both lenses are exactly aligned 
along the optical axis. As a second lens is introduced on the op- 
tic axis, the existing Einstein ring generally bifurcates. If the more 
distant of the two lenses is also luminous, then it produces a third 
Einstein ring. In the case of two singular isothermal sphere lenses, 
then the angular radii of the Einstein rings are given by equations 
( |26t , \21i and ( 1281 ). In the case of two point mass lenses, three Ein- 
stein rings arise whose angular radii are given by ( 1371 ). ( |38| ) and 
( 1391 ). These expressions are our second result, the main point being 
that a single source can, in fact, give rise to more than one Einstein 
ring. 

We also discussed briefly the image configurations for these 
two models if axisymmetry is broken. In the case of the two sin- 
gular isothermal lenses, up to two cut circles and four images can 
arise. In the case of the two point lenses, we find that there are 
always four images. In both instances, the images turn out to pos- 
sess lensing invariants - for example, the magnification invariants 
([33} and ( 142 b . This has not been established before for multi-plane 
lensing and is our third result. 

Finally, we make some critical remarks about possible exten- 
sions of this work. After noting the existence of multiple Einstein 
rings due to a single source, the next obvious question is to ask how 
many there can be. But counting Einstein rings is clearly more diffi- 
cult than counting discrete images. This is because discrete images 
are non-degenerate stationary points of the time delay surface, so 
theorems of Morse theory can be used, whereas Einstein rings are 
degenerate stationary curves and Morse theory does not apply. Al- 
ternatively, the problem of counting Einstein rings can be set up as 
a one-dimensional fixed point problem according to equation ( |17t 
within the framework of intersection theory in algebraic geometry 
and topology. However, there is a degeneracy problem here, too, 
because a non-transverse intersection occurs if the Einstein ring is 
degenerate in the sens e that, in one lens plane, k = 1 at this radius 
dSchneider et al. 1999I . p. 234). Given the large number of free pa- 
rameters in the case of multiple lens planes, it will be difficult to 
ensure the absence of degenerate Einstein rings. 



REFERENCES 

Belokurov V. et al., 2008, MNRAS, submitted 
('arXiv:0806.4188':astro-ph] ) 

Bronshtein I. N., Semendyayev K. A., 1985, Handbook of Math- 
ematics, Verlag Harri Deutsch, Frankfurt am Main 
Eigenbrod A., Courbin F, Dye S., Meylan G., Sluse D., Vuissoz 
C, Magain P., 2006, A&A, 451, 747 
Einstein A., 1936, Sci, 84, 506 
Erdl H., Schneider P., 1993, A&A, 268, 453 
Gavazzi R., Treu T., Koopmans L. V. E., Bolton A. S., Moustakas 
L. A., Buries S., Marshall R I., 2008, Apl, 677, 1046 
Hunter C, Evans N.W., 2001, Apl, 554, 1227 
Kochanek C. S., Apostolakis J., 1988, MNRAS 235, 1073 
Nemiroff R. J., 1993, Am. I. Phys., 61, 619 
Fetters A. O., 1995, 1. Math. Phys., 36, 4263 
Schneider P., Ehlers J., Falco E. E., 1999, Gravitational Lenses, 
Springer- Verlag, Berlin 

Werner M. C., 2007, 1. Math. Phys., 48, 052501 
Witt H. J., Mao S., 1995, ApJ, 447, L105 



ACKNOWLEDGMENTS 

MCW gratefully acknowledges financial support by the STFC of 
the United Kingdom. JA thanks for the hospitality during the vis- 
its to Cambridge. JA also acknowledges that the Dark Cosmol- 
ogy Centre is funded by the Danish National Research Foundation 
(Danmarks Grundforskningsfond) . 



© 0000 RAS, MNRAS 000,[T]l7] 



